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We consider fermionic chains where the two halves are either metals with different band- 
widths or a metal and an insulator. Both are coupled together by a special bond. We 
study the ground-state entanglement entropy between the two pieces, its dependence on 
the parameters and its asymptotic form. We also discuss the features of the entanglement 
Hamiltonians in both subsystems and the evolution of the entanglement entropy after joining 
the two parts of the system. 


I. INTRODUCTION 


The entanglement properties of free-fermion systems have been the topic of many studies and 
various different cases have been investigated, see e.g. [l|. In one dimension, these comprise 
homogeneous critical chains where the ground-state entanglement entropy S varies logarithmically 
with the size L of the subsystem, and non-critical ones where it ^proaches a constant for large 
L. Single defects, both at interfaces 


inhomogeneous systems with random 
couplings or random site energies 


2-lC| and in the interior ll-ldl were studied, as well as 

ui n n□ 

l44l7j|. aperiodic [1^ and exponentially decaying [19|, l20(] 


2ll-l23l|. Finally, systems in external potentials which produce a 


varying density and surface regions have been investigated 


24 


30(]. 


In the present work, we look at yet another situation, namely at systems formed from two 
pieces with different properties. Specifically, we study chains composed either of two critical parts, 
or a critical and a non-critial one. Such systems have been considered previously in the context 
of conformal invariance 3ll-l34l|. In our case, they are realized in the form of undimerized or 


dimerized tight-binding models coupled by a special bond. Physically, this corresponds to either 
two metals with different bandwidths, or a metal and an insulator, and we will use this terminology 
in the following. In both cases one has two types of single-particle eigenfunctions: those essentially 
confined to one of the subsystems, and other ones extending through the whole system but having 
different wavelengths in the two parts. This is the same situation as for a potential step in quantum 
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mechanics. The occupied single-particle eigenfunctions determine the entanglement, and we study 
it between the two different pieces of the composite system in its ground state. 

For the metal-metal system at half filling, the asymptotic result is very simple. It turns out that 
only the interface bond matters and one comes back to the defect problem solved previously in [4 ]. 
Thus S varies logarithmically and the coefficient Ceff is determined by the transmission through the 
interface. The subleading terms, on the other hand, depend on the asymmetry but one can take 
this into account by a rescaling of the length and find links to conformal formulae based on the 
nature of the extended states. Away from half filling, the entanglement is small as long as only the 
localized states are occupied. It increases, when the extended states come into play, but there are 
strong variations with the filling which depend on the size and are also connected with the ratio 
of the bandwidths, which brings an additional length scale into the problem. 

For the metal-insulator system at half filling, the entanglement lies systematically between that 
of the two pure systems. At large sizes, the insulator with its band gap dominates and limits the 
increase of S. For sizes smaller than the correlation length, however, there are no states in the gap 
and one observes a logarithmic increase of the entropy. The system is also a good case to compare 
the entanglement Hamiltonians on both sides. They have the same spectra, but their single-particle 
eigenfunctions and their explicit forms as hopping models are quite different. Basically, one finds 
that the features are similar to those of the pure systems on the corresponding side of the chain. 
For the eigenfunctions this means a decw from the interface into the interior which is slow in the 
metal and exponential in the insulator |ll. l35|. 

For the metal-metal system we also study the behaviour of the entanglement entropy after 
a local quench in which the two pieces in their ground states are put together. Here the two 
Fermi velocities can be seen directly in the time structure and the result can be interpreted in the 
well-known picture of two emitted particles 3G|. 


The paper is organized as follows. In section 2 we describe the set-up und give the basic formulae. 
In section 3 we investigate the metal-metal case by forming and diagonalizing numerically the 
correlation matrix. We show the entanglement entropy and discuss its scaling behaviour and filling 
dependence together with some entanglement spectra. In section 4 we consider the metal-insulator 
case at half filling and compare all relevant quantities with those of the pure systems. In section 5 
we present the time evolution of S after connecting two metallic systems and Section 6 contains a 
summary. Finally, some analytical results for the metal-metal system are given in an appendix. 
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II. SETTING AND BASIC FORMULAE 


We consider open chains of free fermions with nearest-neighbour hopping and 2L sites. The 
hopping is different in the left and right half and also between both parts. The Hamiltonian is 


L-l 


H — in (chcn+1 + 


n=—L+1 


For the metal-metal system 


( 1 ) 


tn = 


tl 

to 

t2 


n < 0 
n = 0 
n > 0 


( 2 ) 


and one has three parameters in the problem. However, by a rescaling of H, one can always 
achieve tit 2 = 1. We will assume this in the following and write the quantities as ti = exp(A), 
t 2 = exp(—A) and to = exp(Ao). Moreover, we will always consider positive A, i.e. ti > t 2 - In 
some places we also use the ratio r = t 2 /ti. 

For vanishing coupling (to = 0), the two parts of the chain have single-particle energies 


OJa = —taCOs{qa), a = 1,2 


( 3 ) 


with momenta qa = 'Km/{L + l),m = 1,2,...L. The corresponding band structure with bands 
between ± ta is shown on the left of Fig. [TJ In the coupled system, one has extended states for 
|a;| < t 2 , while outside this region they are confined essentially to the left half-chain. Moreover, 
one can have two states localized at the interface if fo is large enough. 

For the metal-insulator system we take 


1 : n < 0 


tn. — ^ 


to : n = 0 


( 4 ) 


1 + (-1)^+M : n > 0 

In the insulator, one thus has alternating hopping, two sites per unit cell and the single-particle 
energies 


C02 


Z\J cos2(g2) -I- (5^ sin^(g2) 


( 5 ) 


where the momenta have to be determined from the boundary condition. Thus there is a gap 
between ±<5. The resulting band structure is shown on the right of Fig. [TJ For an open chain, 
there are also states localized at the boundary with exponentially small a; 2 , if the outermost bond 
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FIG. 1: Schematic band structure of the uncoupled composite systems. Left: Metal-Metal, Right: Metal- 
Insulator. 

is a weak one. The coupled system has extended states for |a;| > d and others, confined essentially 
to the metal, inside the gap. In this sense, it is still critical. 

It is simple to set up the eigenvalue problem for the full systems, and some details are given in 
the Appendix. However, the matching conditions in the center are difficult to handle in general. 
Therefore we determine the single-particle eigenvalues ojq and the corresponding (real) eigenfunc¬ 
tions 4>g(n) numerically. The correlation matrix Cmn = (cmCn) is then 

Cmn = ^ 4>q(m)4'q(n) (6) 

qeF 

where the sum extends over the states q in the Fermi sea. Restricting Cmn to the chosen subsystem 
(either the left or the right half-chain), its eigenvalues Cfc give the single-particle eigenvalues = 
ln[(l — Ck)/Ck] of the free-fermion entanglement Hamiltonian T-L in 

p = ^exp{-n) (7) 

where p is the reduced density matrix [^. From them, the entanglement entropy S follows as 

^ = Eln(l+e--) + E^ (8) 

k k 

In addition to the eigenvalues Cfc resp. we also determine the corresponding eigenfunctions pk 
and construct H in section 4. 


III. METAL-METAL SYSTEM 

We first consider half-filled systems where the Fermi level is in the middle of the bands. The 
resulting entanglement entropies for Iq = 1 and three values of the parameter A are shown in Fig. 
[21 both for even and for odd L. One sees that the typical increase with the size, known for the 




















5 


homogeneous case A = 0, persists in the composite systems. The plot against InL shows that 
also the asymptotic law l/61nL + /c is unchanged. However, the value of k becomes smaller and 
the finite-size effects, in particular the even-odd alternation, increase dramatically as A becomes 
larger. For A = 3 (not shown), the values of S are initially close to zero and to In2, respectively. 
In this case, the ratio of the bandwidths t 2 /ti ~ 1/400 is already very small, the states localized on 
the left drop rapidly on the right, while the extended states have very small amplitudes on the left 
(see the Appendix). For even L, this leads to a small entanglement which only increases as more 
states come into play at larger L. The value In 2 for odd L is a consequence of the particle-hole 
symmetry of the problem which, in this case, forces one of the Sk to vanish. 




FIG. 2: Entanglement entropy S for Aq = 0 and three values of A. Left: as function of L. Right: as 
function of InL. The upper (lower) curves correspond to odd (even) values of L. 

The function S{L,A) has an intriguing universal behaviour. If one defines a length L* = L/a 
and chooses the scale factor a properly, one can achieve S{L*, A) = S{L, 0), i.e. all curves collapse 
on the one for the homogeneous system. Written differently, one has 

5(L,A) = 5(aL,0) (9) 

In the asymptotic region, where one has straight lines in a logarithmic plot, this feature is not sur¬ 
prising: A rescaling can always generate the necessary shift to make the curves coincide. However, 
the phenomenon is not restricted to this region. As shown in Fig. [3] on the left, it also holds for 
small sizes where there is still curvature in the graph. The variation of a with A is given on the 
right hand side of Fig. [3l One sees that the results determined from different (even) L agree very 
well and that a becomes rapidly smaller as A increases. The curve can be described approximately 
by 1/ch^ A. Note that the relation ([9]) can only be applied to values of L such that aL is larger 


than 1. 
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FIG. 3: Scaling behaviour of the entanglement entropy. Left: S' as a function of the effective length aL for 
four values of A (logarithmic scale). Right: Variation oi L/L* = a with A. The dotted line is the conformal 
factor in cni). 


At this point, an observation from the appendix is helpful. Namely, near the Fermi level the 
wavefunctions in the right subsystem are the same as for a homogeneous system of total length 
L(1 + r). Thus, ignoring the states outside the narrow band, the correlation matrix on the right is 
the same as for an unsymmetrically divided homogenous chain and one can invoke the conformal 


result for the entropy 


38l | extended in [39| to include 1/L corrections. This gives 


1 (- 1 )^ 

Spc = - Inz-h const., 

6 z 


4 TT^ 

z = —L[l + r) sm —---. 

vr ^ ’ L{l + r) 


( 10 ) 


with i = L. The r-dependent factor in z provides a similar rescaling of L as a, and a corresponding 
shift of the curves. It turns out, however, that this is not enough to make them coincide. Rather, 
one has a residual difference S — Spc = Sj{A), which one can attribute to the interface, and which 
rises smoothly from 0 to about 0.13 as A increses from 0 to 3. It has its origin in the neglected 
exponentially decaying states below the narrow band. 

This can be seen clearly, if one considers a partition of the system with i < L sites on the right 
and calculates the entanglement between this subsystem and the remainder. As demonstrated in 
Fig. H the expression (|10l) then fits the data very well as soon as one moves away from the interface 
by a few lattice sites. One can also consider a division located in the left part of the system where 
the extended wave functions are the same as in a homogeneous system with total length L(l + l/r), 
but in this case the conformal formula does not work well because the additional states are more 
important. These results also show that the symmetrical division is actually a somewhat special 


case. 


We now turn to the effect of to- If one varies the central bond in a homogeneous chain (A = 0) 
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FIG. 4: Entanglement entropy between the outermost I sites of the right subsystem and the remainder for 
L = 100 and three values of A. The inset shows the deviation from the formula oni). 


the asymptotic behaviour of S is 


5 = — InL + A: 
6 


( 11 ) 


where Cefr depends on to and is given by an explicit, tho^h lengthy formula in which only the 
transmission coefficient at the Fermi energy T = enters [^. But the calculation in the Appendix 
shows that T is independent of A and given by 

1 


T = 


ch^ Ar 


( 12 ) 


Therefore one expects the result of the homogeneous problem also in the composite systems. This 
is, in fact, the case and demonstrated in Fig. [S] on the left, where results obtained by htting the 
data between L = 100 and L = 400 to ifTT]) plus a 1/L term are collected. For A = 2, the average 
between even and odd sizes was taken, while for A = 1 this was unnecessary. This verifies the 
Fermi-edge character of Ceff also in the present case. 

The constant k, on the other hand, depends on A, as was found above already for the special 
case Iq = 1. It is shown in Fig. [5] on the right. The values for odd L are always larger than those 
for even L with the maximal difference of In 2 appearing for to —)• 0. This is connected with the 
half-integer particle number in each subsystem for odd L which leads to a kind of singlet state 
even in the decoupling limit. The minimum in k for even L is also present for A = 0 and was also 
found for a segment in an infinite chain Qj. 
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FIG. 5: Effective central charge CeS and constant k in ([IT 
line A = 0 in the left figure is the theoretical result, see |4| 


as functions of to for different values of A. The 


Finally, we consider the dependence of S on the filling v = N/2L of the system, where N is the 
number of particles. It is shown in Fig. [6] for a fixed length L = 200 and four values of A. The 
smooth curve (red) is the result for the homogeneous system and given, up to the small even-odd 
oscillations 39-^ and a constant, by l/61n(singi7’), where qp = [^]. Turning on A, regions 


of very small entanglement develop for small and large filling, where only states localized on the 
left are occupied or remain empty. At the same time, the oscillations of S in the central region 
become slower and slower. 
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FIG. 6: Entanglement entropy S' as a function of the filling ly = N/2L for L = 200 and different A. 

The origin of these oscillations, at the level of the eigenvalues can be seen in Fig. [71 There, 
the low-lying are shown for all fillings between 0.2 and 0.8. For the case A = 2, shown on the 
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right, the variation of the with v is rather slow and the curves cross zero only at two fillings. At 
these points, the crossing eigenvalue gives the maximal contribution to 5, namely In 2, and they 


coincide with the locations of the maxima in Fig. 
much faster and one has altogether 16 crossings, 
center of Fig. [6l 



V 


El For A = 1, shown on the left, the variation is 
This is again the number of maxima of S in the 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 


V 


FIG. 7: Variation of the low-lying single-particle eigenvalues with the filling v = N/2L for L = 200. Left: 
A = 1. Right: A = 2. The colours are the same as in Fig. El 


These structures can again be understood from the relation to a smaller, unsymmetrically 
divided homogeneous system. In such a case, if L' is the size of the smaller subsystem, one has 
L' crossings of the eigenvalues and the effect on S is given by a factor sm{qF(2L' + 1))/sin^i? 
replacing the alternating sign of the Xjz term in (|inii . This leads to slow oscillations, although less 
pronounced than the observed ones. However, L' = rL gives L' = 3.6 for A = 2 and V = 27 for 
A = 1 and thus not the correct numbers (2 and 16), even if one takes the nearest integer. But 
one can argue that the length rL only refers to the states in the center of the band, while near the 
edges sinA :2 appears in (1261) which is smaller than k 2 by a factor of 2/7r. Working with L' = 2rLl'K 
gives values 2.3 and 17.2, respectively, which are quite close to the numerical findings. One also 
arrives at U by counting the number of levels which the left subsystem contributes to the total 
number of band states. To accommodate them, one needs just this number of additional sites. 

One should mention that one can see basically the same structures also in the energy u) as 
tiny deviations from the otherwise smooth level spacing of the extended states. In this case, the 
interpretation is that L' also gives the number of periods of the slowly varying tangent in (|26p as 
k 2 sweeps through the band, and thus the number of possible peculiarities in the allowed momenta. 
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IV. METAL-INSULATOR SYSTEM 

In the following, we always consider half-filled systems. To avoid boundary states in the dimer¬ 
ized subsystem, we work with even L and choose strong bonds 1 + <5 at the borders. The entan¬ 
glement entropies for two relatively small dimerizations and to = 1 are shown in Fig. [SJ together 
with those for the pure systems. One sees that the mixed system always has an intermediate 
entanglement, which is very plausible. Replacing one half of the metal by an insulator reduces the 
metallic result, while replacing one half of the insulator by a metal increases the insulator result. 
For 5 = 0.001, shown on the left, all three curves look similar and there is no sign of a saturation. 




FIG. 8: Entanglement entropies for metal-metal, metal-insulator and insulator-insulator systems as a func¬ 
tion of L for two values of the dimerization parameter S. 

while for 6 = 0.01 a saturation is clearly visible for the lower ones. This can be interpreted in 
terms of the correlation length ^ = 1/6 in the insulator. It is ^ = 1000 in the first case so that the 
non-criticality does not fully show up. In the second case, where ^ = 100, it does and the crossover 
takes place roughly around this value. 

In the region L one can fit the curves for the composite system to a form a In L+b and finds 
values of a quite close to 1/6. This can also be done for to < 1 and gives curves for c^s which lie 
somewhat below the metallic one shown in Fig. [5j but approach it as 5 goes to zero. In the opposite 
limit, L 3> the entropy saturates and for to = 1 one finds the asymptotic law S ~ 1/6 In ^-|-const, 
as for the pure insulator but with a larger constant. To obtain the coefficient 1/6 accurately from 
not too large one has to include sub leading corrections of the form (dln^ -|- e)/^ which one 
knows in the pure case from the connection to a transverse Ising chain sketched below. Again, the 
analysis can be extended to to < Ij and in this case one finds c^s to an accuracy of 3-4 digits. 

On the level of the wavefunctions entering the correlation matrix, the (pseudo-) critical be- 
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haviour can be understood in the following way. For L < ^ one has \u\ ^ 1/L > 5 and there are 
no states in the metal with energy in the band gap of the insulator and none in the insulator too 
close to the band edge. Thus one only has extended wavefunctions which are similar to those in 
the pure metal. This can be seen analytically or from the numerics. 



FIG. 9: Eigenvalues Ek for metal-metal, metal-insulator and insulator-insulator systems with L = 50 and 
(5 = 0.1 for the insulator. 


It is also instructive to look at the £k spectra. In Fig. [9] an example is shown for the case 
(5 = 0.1 and L = 50. The lowest curve with the slight bend is the well- known result for the metal, 


while the hh 


rest one is for the insulator and strictly linear as for the transverse Ising chain or the 


XY model [ij. In fact, one can find a relation between the for the dimerized hopping model 


and the transverse Ising chain using the method in 
coupling A = (1 — (5)/(l -|- (5) and the £k are 


44l |. The latter model then has field h = 1 and 


ek = {2k + l)£, k = 0,±l,±2,..., e = n I{X')/I{X) 


(13) 


where /(A) denotes the complete elliptic integral of the first kind and X' = \/l — A^, see also |45( |. 
The MI result lies in between, with a small break after the first eigenvalue, which becomes more 
apparent for larger values of <5, when the curve becomes steeper and more linear. 

For the spectrum, it does not matter in which of the two subsystems one considers the correla¬ 
tion matrix. Formally, this follows from the property = C for the full system [^ . However, the 
eigenvectors of the two reduced matrices, and therefore the corresponding entanglement Hamilto¬ 
nians Ti, are quite different. This is illustrated in Fig. [TOlfor the lowest \ek\ and 6 = 0.1. We have 
plotted ^k{n) for +£k on the left and for —on the right, since these two are directly related. 
Namely, if one partitions C into four blocks according to the location of the sites and is an 
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eigenvector of for part 1 with eigenvalue C; then 

is an eigenvector of for part 2 with eigenvalue 1 — C- 

One sees the typical features of the pure systems on the two sides, namely a slow power-law 
decay in the metal (left) and a fast one within a distance ^ = 10 in the insulator (right). The pure 
cases are shown for comparison, and one notes that the differences in the composite system are 
relatively small, lying more in the amplitudes than in the general behaviour. 
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FIG. 10: Eigenvectors of the Ml correlation matrix for the smallest |efc| in the metallic (left) and the 
insulating (right) subsystem for L = 50 together with the results for the homogeneous cases (MM and 11). 


With the eigenvectors, one can explicitly determine the entanglement Hamiltonians %. They 
have again the form of hopping models with basically only nearest-neighbour hopping as in H, but 
the amplitudes tn increase from the center towards the boundaries. 

In the pure metal, this increase is monotonous and the in vary as n(2L — n) [l|, while in the 
insulator one has an additional alternation which is coupled to the bond alternation in H. Due to 
the form of the eigenvectors, the result in the composite system is again close to that of the pure 
system on the corresponding side. This is shown in Fig. [TT] for a system with a total of 2L = 16 
sites, which is the largest size attainable before numerical errors in the large £k set in. 
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FIG. 11: Nearest-neighbour hopping in the entanglement Hamiltonians for the left and right subsystem as 
a function of the position together with the results for the pure systems for L = 8 and (5 = 0.1. 


V. EVOLUTION AFTER A QUENCH 


In this section, we return to the metal-metal systems and study how the entanglement evolves 


after one joins the two initially disconnected parts. This type of local quench has been considered 


repeatedly in the past 


| n |,|47 


-j59l|. For the homogeneous case, one finds oscillations of S{t) which 


one could call “entanglement bursts”, see e.g. 


5ll |. In our case, the situation is more complex, 


because of the interface and the two Fermi velocities. 

The two pieces are coupled at time t = 0 by switching on to and the evolution of C{t) is 
calculated via the Heisenberg operators Cn{t), ch{t). The resulting S{t) is shown in Fig. [12] for two 
different values of A and of to- The system on the left with A = 0.5 has roughly the band widths 
of Fig. [ll while on the right the asymmetry is larger. 




FIG. 12: Evolution of the entanglement entropy S after connecting two metallic systems with L = 100. 
Left: A = 0.5. Right: A = 1. The times indicated are the theoretical values, see the text. 
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For Aq = 0 (upper curves) one has the following general features: An initial burst ending with 
a decrease at a time ri followed by a plateau, another burst at time T 2 which is a kind of mirror 
image of the first ending at time ti + T 2 and then a rough repetition. For A = 1, the plateau is 
longer and shows six additional structures, whereas for A = 0.5 only a single one is visible. The 
lower curves have roughly the same pattern, but apart from the first maximum the structures are 
more washed out. 


These features can be interpreted in the well-known picture due to Calabrese and Cardy 


361, in 


which a pair of particles is emitted from the junction at time t = 0, travelling in opposite directions 
and spreading the entanglement. In the present case, they have the velocities vi = ti and V 2 = t 2 
and the corresponding space-time diagram is given in Fig. fTHl 



FIG. 13: Trajectories of the particles in the Calabrese-Cardy picture. 


After time ri = 2L/vi, the left particle returns to the center, and if it simply moved into the 
right subsystem, the entanglement would continue to drop. This happens in a somewhat related 
situation where the velocities are the same, but the size of the right subsystem is larger Here, 


however, there is a probability for reflection at the interface, such that the entanglement rises again 
and another small “burst” follows. For A = 1, where vi = 7Av2, the left particle can make seven 
round trips before the right particle retuns to the center at time T 2 = 2L/u2 = Ti/r. These are very 
clearly visible in the Figure. There are, however, some remarks to make. According to (jl2p there is 
no reflection right at the Fermi level for Aq = 0. The effect can therefore only come from somewhat 
slower particles away from it. In fact, the numerical values for ri and T 2 are somewhat larger than 
the theoretical ones, which supports this argument. Also, one would expect the reflection effect to 
become weaker with each cycle, which is only barely the case. On the other hand, for Aq = — I 
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one has R = 0.58 at the Fermi level and the effect should be stronger, which is also not the case. 
Nevertheless, the picture gives a good overall description. 

One could presume that at least the first burst can be described by the conformal result 


m 


Sh{vt) = 2 In 


2L TTUt 
— sin —— 
vr 2L 


+ const. 

valid for a homogeneous system with particle velocity v. Thus we tried the Ansatz 


(15) 


S{t) = [Sh{vit) + Sh{v 2 t )]/2 + const. (16) 

which, on the scale of Fig. 12, fits the data for Aq = 0 very well. However, a closer look at the final 
decrease shows that it comes too early since, as mentioned, the value of ri is too small. Finally, 
we note that the first-burst results for a homogeneous but unequally divided system with the same 
round-trip times ri = 2Li and T 2 = 2 L 2 lie above ours, and the difference increases with A. Thus 
the two problems are not trivially connected. 


VI. SUMMARY 

We have considered the entanglement in fermionic chains composed of two different halves. In 
solid-state terminology, they were metals and insulators, while in statistical-physics terms they 
corresponded to critical and non-critical systems. The metal-metal system at half filling showed 
the same logarithmic behaviour of S' as a homogeneous chain with a defect. In the subleading 
terms, the asymmetry showed up, but the entropy had a remarkable scaling property. Its filling 
dependence, finally, showed oscillatory behaviour again linked to the asymmetry. Both features 
could be understood from the nature of the extended eigenfunctions which are the same as in 
homogeneous systems with different lengths. 

The difference in the Fermi velocities showed up directly in the quench experiment where two 
half-filled metals were joined and the entanglement was monitored. In that case, the value of 
r determined, how many cycles one sees in S{t), before a certain return to the initial situation 
takes place. This pattern would occur in unsymmetrically divided homogeneous systems only with 
additional defects. 

The metal-insulator system was somewhat simpler. Its entanglement properties were seen to 
lie always between the two pure systems and as a function of the size, one has a crossover from 
an initial critical behaviour with a logarithmic increase of the entanglement to a saturation typical 
for a non-critical system with finite correlation length. It also provided an instructive example 
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for entanglement Hamiltonians which are quite different for the two subsystems although they 
have the same spectra, as they must. One should note that the dimerized hopping model used for 
the ^ulator, and investigated already in plays a central role in the theory of polyacetylene 


60, ISlI]- In that sense, we studied a particular polymer system. 

We considered the geometry with open ends, because then one has only a single interface. 
However, one could equally well look at rings. For composite transverse Ising models this was done 
first in |3ll . l32l | but without particular interface bonds and with a view on the spectra. Obtaining 
the entanglement entropy is more complicated, and we could only give an analytical result for the 
asymptotic behaviour of the metal-metal system. 

While we focussed on the entanglement entropy, it is known that the particle-number fluctua¬ 
tions in the subsystems behave similarly 


62l |. Thus, in a homogeneous metal, the prefactor of the 
InL term in the fluctuations is given by T/27r^ with the transmission coefficient T Q] and varies 
qualitatively like The same result is found in the composite systems, and this could offer a 
way to access the entanglement experimentally. 
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Appendix: Some formulae for the metal-metal system 


Let (j){n) and </>(n) denote the wave function at site n on the left and right, respectively. The bulk 
solutions on the two sides are waves with momenta qi and q 2 and energy io = —ti cos qi = —t 2 cos q 2 - 
At the interface, the equations are 

to ^{0) = t2 4>{0) , h = to (17) 


Eor a wave coming in from the left and being partially reflected and transmitted, the amplitudes 
are 


4>{n) = Ai exp (iqin) + Bi exp {—iqin) 
(j){n) = A 2 exp {iq 2 n) 


( 18 ) 
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Inserting this into the relations (1171) . one finds for the reflection coefficient 


R = 


Bi 


ch 2i^ — cos (gi — ( 72 ) 


(19) 


ch 21 / — cos {qi + ( 72 ) 

where the quantity u is defined as exp {2i/) = tQ/tit 2 - In our parametrization, where tit 2 = 1, one 
has 1 / = Aq. The transmission coefficient follows from T = 1 — R and can be written 


T = 


sm qi sm q 2 


( 20 ) 


sh^ u + sin^ ((gi + 92)72) 

If the left and right subsystems are identical, 91 = 92 = Q, and one finds the result for a bond defect 
in an otherwise homogeneous chain. The same holds in a general system if one is in the middle of 
the bands. Then w = 0, 91 = 92 = vr/2 and one has the result 

1 


T = 

ch^ 1 / 

Due to the factor sin 92 in (j20p . the transmission vanishes at the edges of the narrow band. 
The extended eigenfunctions of H have the form 

(7(n) = Ais\nqi{n + L) 

(p{n) = A 2 sinq 2 {n — L — 1) 


( 21 ) 


and the conditions (fT71) lead to 


sin9i(L + 1) sin92(L + 1) 


0 

iii2 


( 22 ) 


(23) 


sin9iL sin92L 

which, together with ti cos 91 = t 2 cos 92 , determines the allowed momenta. Setting q^ = vr/2 + k, 
and assuming L even, it takes the form 


cosA;i(L + l) cos/c 2 (T + 1) 


sinkiL sin/c 2 T 


0 

hh 


(24) 


The relation between the reads 


sin fci = r sin ^ 2 , v’= ^ 2/^1 

and for small r simplifies to ki = r sin ^2 • Inserting this into 
numerator then gives 

t2 

cot (/C 2 T) = —|;a,n {rL sin k 2 ) 
tlt2 


(25) 


and setting L + 1 ~ L in the 


(26) 


which contains only the momentum k 2 - The solutions can be discussed graphically, but one sees 
directly that besides L a second length L' = rL appears. If it is large enough, the tangent 
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completes several cycles as k 2 varies, and this can lead to additional features in the momenta or in 
the amplitude ratio which is, again with L + 1 ~ L, 

M _ _h cos(rLsin^ 2 ) 

Ai to smk2L 

If /c 2 is small, ki = rk 2 for all r < 1 and (1261) for t^/tit 2 = 1 can be written in the two forms 

tan (^ 2 ?*-^) tan (/c 2 -h) = 1 , tan (/siL) tan (/ciL/r) = 1 (28) 

Both equations are the same as for a homogeneous chain built up from two pieces with different 
lengths and can therefore be solved explicitly. For k 2 , the homogeneous system has length L + L', 
while for ki it has length L + L", where L" = L/r. Thus the small momenta are equidistant, i.e. 

A :2 = ^^^^(2n + 1), n = 0,±l,±2,... (29) 

Using ([25]) in (f27)l . one finds {A 2 /A 1 Y = t\ = l/r, i.e. the amplitudes are larger on the right than 
on the left. Including normalization to leading order, this gives prefactors Af = 2/(L + L") and 
A 2 = 2l{L + L') which correspond exactly to these chain lengths. 
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